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[JLh ! Abstract 

■ Nuclear-medium effects in the weak structure functions F2{x,Q'^) and F3{x,Q'^) in the charged 

current neutrino and antineutrino induced deep inelastic reactions in 208p|-, j^^ve been studied. The 

calculations have been performed in a theoretical model using relativistic nuclear spectral functions 

^ ■ which incorporate Fermi motion, binding and nucleon correlations. We also consider the pion and 

a' 

rho meson cloud contributions calculated from a microscopic model for meson-nucleus self-energies. 
. Using these structure functions, the results for the differential cross section have been obtained and 

I compared with the CERN Hybrid Oscillation Research apparatus (CHORUS) data. The results 

^ ; for the ratios -^^^ 208F^'= ' wsh' ' 208F-'^ ' ^'^'^ '^F^ (i=2,3) have also been obtained and a few 

have been compared with some of the phenomenological fits. 
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I. INTRODUCTION 



In recent years the need for a better understanding of nucleon dynamics in the nuclear 
medium has been emphasized upon in the weak interaction induced processes in order to 
precisely estimate the lepton event rates in the Monte Carlo generators which are being 
used in the analysis of present generation of neutrino oscillation experiments [l,l2|- Most of 
these experiments are using neutrino and antineutrino beams in the few GeV energy region 
and the targets are nuclei like ^^C, ^^Fe, ^°^Pb, and so on. In the few GeV energy region 
the contribution to the cross section comes from the quasielastic, inelastic as well as the 
deep inelastic processes. In the deep inelastic region, both experimentally as well as theo- 
retically, limited efforts have been made to understand the medium effects in the structure 
functions. Substantial number of experiments have been performed using charged lepton 
beams to study nuclear modifications of the electromagnetic structure function F^^{x, Q^) 
with several nuclear targets while few experiments have been performed using neutrino and 

n n 

antineutrino beams on some nuclear targets like carbon, neon, iron and lead m-M- In the 
case of experimental measurements performed for F^^^"'''{x, Q"^) the error bars are large and 
need better precision. In general, nuclear modifications in the weak structure functions 
F2^*^"'^(x, Q^) may be different from the nuclear modification in F,^'^^{x,Q'^). Particularly 
neutrino data are important in the determination of valence quark distributions in the nu- 
cleon. This is due to the fact that the parity-violating (^p^'^"-^^ structure function directly 
probes into the valence distributions. Precise determination of parton distribution functions 
(PDFs) is also necessary for the new physics at the colliders. But for the determination of 
PDFs, nuclear medium effects should be properly accounted for. Nuclear-medium effects 
may be responsible for the anomaly observed in the NuTeV experiment for the weak mixing 
angle sin^ 9w js]. NuTeV Collaboration js] has performed measurements of neutrino and 
antineutrino induced deep inelastic scattering cross sections using iron target. The dif- 
ferential cross sections are then used to extract weak structure functions F^'^°'^{x,Q'^) and 
pWeak^^^Q2y CHORUS Collaboration jo, 7| has performed high-statistics measurement of 
the differential (anti)neutrino cross sections using mainly lead target at various neutrino and 
antineutrino energies as a function of Bjorken scaling variables x and y. Analysis is being 
made at the Neutrino Oscillation MAgnetic Detector (NOMAD) jsl for the weak structure 
functions and cross section measurements with carbon target using neutrino beam. The 



Neutrino Scattering On Glass (NuSOnG) experiment 10| has been proposed at Fermilab to 
study the structure functions in the deep inelastic region using neutrino scattering on carbon 
target. Main INjector Experiment v-A (MINERz/A) ll| is taking data using neutrinos from 
the Neutrinos at the Main Injector (NuMI) facility, measuring neutrino cross sections in the 
energy region of 1-20 GeV using various nuclear targets and the aim is to study the strong 
dynamics of the nucleon and nucleus that affect these interactions as well as the parton dis- 
tribution functions. The Oscillation Project with Emulsion-t Racking Apparatus (OPERA) 
experiment jjjl is a long baseline experiment using lead emulsion target and the main pur- 
pose of the experiment is the observation of to Pr oscillations in the direct appearance 
mode. The charged current deep inelastic scattering process is found to be dominant in this 
experiment with a fraction exceeding 90%. Therefore, several experimental activities are 
going on to study neutrino as well as antineutrino event rates in the deep inelastic region. 
On the theoretical side, for the weak interaction induced processes in the deep inelas- 



tic region, t 
authors 13- 



le dynamical origin of the nuclear medium effects has been studied by a few 



15| . In some other theoretical analysis, nuclear medium effects have been phe- 



nomenologically described in terms of a few parameters which are determined by fitting 



the experimental data o: 
various nuclear targets 



c harg ed leptons and (anti)neutrino deep inelastic scattering from 



16 



23(1 . The various phenomenological studies differ in the choice 



of the data sets(Lepton-Nucleus-|-Drell Yan data, Lepton-Nucleus+Drell Yan+z/(z/)-Nucleus 
data, and z/(z/)-Nucleus data), experimental cuts in their analysis, parametrization of the 
parton distributions, etc. 

In the present work, we study nuclear medium effects on the structure functions F2{x,Q'^) 
and F^{x,Q'^) in lead, treating it as a nonisoscalar nuclear target. This study has been per- 



formed using a relativistic nuclear spectral function [2J] to describe the momentum distribu- 
tion of nucleons in the nucleus within a field-theoretical approach where nucleon propagators 
are written in terms of this spectral function, and nuclear many-body theory is used to cal- 
culate it for an interacting Fermi sea in nuclear matter. A local- density approximation is 
then applied to translate these results to finite nuclei Q, [ ' 



26|. 



We have assumed the Callan-Gross relationship for nuclear structure functions Fi^{x^ Q'^) 
and F2'^{x, Q'^). The contributions of the pion and rho meson clouds are taken into account 
in a many-body field-theoretical approach which is based on Refs. Q, l^^] . We have taken 
into account target mass correction (TMC) following Ref. 28| which has a significant effect 



at low Q^, moderate, and high Bjorken x. To take into account the shadowing effect, 
which is important at low and low x, and modulates the contribution of p ion and rho 



cloud contributions, we have followed the works of Kulagin and Petti |l3l . |29| . Earlier we 
have applied the present formalism to study nuclear effects in the electromagnetic structure 
function F,^^\x,Q'^) in nuclei in the charged lepton-nucleus deep inelastic scattering 25| 
as to the study of weak structure functions Fj(x, Q^) (i=2,3) in carbon and iron 



as wef 
nuclei 



15| using (anti)neutrino-nucleus scattering. In the case of electromagnetic interaction 



we have found that the results of our calculations are in agreement with recent results as 
obtained from the Thomas Jefferson National Accelerator Facility (JLab) 30| and also with 
some of the earlier experiments performed using heavier nuclear targets. While in the case of 
neutrino and antineutrino induced reactions, the results in carbon and iron were compared 
with the experimental results of NuTeV js] and the CERN Dortmund Heidelberg Saclay 
Warsaw (CDHSW) fs^ data and found to be in agreement. In this paper, we have studied 
nuclear-medium effects in the neutrino and antineutrino induced deep inelastic scattering 
on lead, and obtained the results for the weak structure functions F^(x,(5^) and F^(x,Q'^). 
Using our results for the structure functions in helium, carbon, oxygen and iron [15], we 
have obtained the results for ^^"^ (i=2,3). These results 

may be quite useful in the analysis of the MINERz/A pUj] experiment. Finally the structure 
functions in lead are used to obtain the differential scattering cross section and the results 
are compared with the experimental results from CHORUS 71]. 

The plan of the paper is as follows. In Sect. [Ill we present in brief the formalism for 
i/(i^)-nucleus scattering, in Sect. IIIII we present and discuss the results of our calculations 
and compare them with the available experimental results. In Sect. IIVI we conclude our 
findings. 



II. FORMALISM 

The differential cross section for charged current neutrino (antineutrino) interaction with 
a nucleus is written as: 

dE' dVL' ^ (27r)2 



k' 

k 



m 



w_ 

m 



w 



al3 T^Mu)A 
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where E' and Q' are the energy and scattering angles {0',<f)') of the outgoing lepton; Gp is 
the Fermi constant of the weak interaction; k' and k are the momenta of the outgoing and 
incoming leptons, respectively; = {k — k'Y is the four momentum transfer square, my/ 
is the mass of the boson W^. and W^j^"^^ are the leptonic and the hadronic tensors 

given by: 

L°J(A;, k') = k'^k'^ + k"^k^ - g"^{k ■k')Ti e"^^^A;^A;^ (2) 



2 e^pp. P'^q'^Wf^^iPAA) + 4l2<lo^^pWf^^^PA.q) 



+ ^ {PAaqp + qaPAp) W^^'^^iPA, q) (3) 

where Pa is the momentum of the nucleus A and W^^{x, Q"^) (i=l-5) nuclear hadronic struc- 
ture functions. Ma is the mass of the nucleus. In the antisymmetric term of the leptonic 
tensor, - sign is for the neutrino and + sign for the antineutrino induced reaction. VF4 and 
W5 are generally omitted in the cross section expression since they are suppressed by a fac- 
tor which vanishes in the limit mi — > 0, being mi the mass of the outgoing charged lepton. 
Finally Wf'i^x^Q'^) (i=l-3) are redefined in terms of the dimensionless structure functions 
F^{x, Q"^) through 

We have also assumed the Callan-Gross relationship for nuclear structure functions Fi^{x) 
and F2'^(x), therefore, we are left with only two independent structure functions viz. F2^{x) 
and F3^(x). 

For the numerical calculations, parton distribution functions for the nucleons have been 
taken from the parametrization of the Coordinated Theoretical-Experimental Project on 
QCD (CTEQ) Collaboration (CTEQ6.6) [si]. The Next-to-Leading-Order (NLO) evolution 



of the deep inelastic structure functions has been taken from the works of Moch et al. 32, 



33| 



In the local-density approximation the nuclear hadronic tensor W^j^^"^ can be written as a 
convolution of the nucleonic hadronic tensor with the hole spectral function. For symmetric 
nuclear matter, this would be: 



where kp{r) is the Fermi momentum for symmetric nuclear matter, depending on the density 
of nucleons in the nucleus, i.e. kpir) = ( ^p(^) ) 

For a non-symmetric nucleus like ^°*Pb, one considers separate distributions of Fermi 
seas for protons and neutrons and the above expression modifies to: 



d^p M 



w:f' = 2 / dv / / dp" sr-'iP". P, k,.,) w:, 



where the factor 2 in front of the integral accounts for the two degrees of freedom of the spin 
of the nucleons. There are two different spectral functions, each one of them normalized 
to the number of protons or neutrons in the nuclear target and are functions of Fermi 
momentum of protons and neutrons respectively which are given by kp^p = (^n'^ppY^^ and 
kF,n = (3vr^p„)^/'^. For the proton and neutron densities in lead, we have used two-parameter 
Fermi density distribution given by 

1 + exp(— ) 

where the parameters cl=6.624 fm, c2=0.549 fm for the protons and cl=6.89 fm and 
c2=0.549 fm for the neutrons have been used js^. We have taken 3-parameter Fermi den- 
sity for Helium, harmonic oscillator densities for carbon and oxygen and 2-parameter Fermi 
density for iron nuclei and the density parameters are taken from Refs. 3^, |35|. The 



oo 



invariant quantities for the deep inelastic scattering (DIS) of neutrinos with nuclei are: 

^2 2 PA-q Pa-Q 

^ =-^' "^=2iV^' "" = ^' '^ = l^k 

where xa is the natural Bjorken variable in the nucleus and xa G [0, 1]; i/a is the inelasticity. 
These two variables are related to the nucleonic ones via: 

Xa = —t'i where x = „ and va = = V (9) 

where x is the Bjorken variable for neutrino-nucleon interaction expressed in the nucleon 
rest frame. We can see that x G [0, A], though for x > 1 the nuclear structure functions are 
negligible. The variable ua varies between the following limits: 

11 

2Ev 2E„ 




0.6 



FIG. 1. (Color online) Ratio R(x,Q 



' 208F^' 



and K{x,Q 



208F, 



using our base equation with 



TMC (dotted line) and the full model (solid line) at LO for = 5 GeV^. The results from Hirai 
et al. [l?! and Eskola et al. 2^ have also been shown. 




therefore, for sufficient high neutrino energy we have < < 1. 



If we express the differential cross section with respect to these variables {xajVa), we 
obtain the following expression in terms of the nuclear structure functions: 



d'a. 



cc 



dxA dyj 
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FIG. 3. (Color online) Ratio R(x,Q^) 



2^_ ^PP 



with full calculation for = 5 GeV^ and 



50 



GeV^ using our base equation with TMC (dashed line) and the full model (solid line) at NLO. 
Here the results for the ratios R(x,(5^)= ^imf' ^imf 208^/-= ^^^^ ^^^^ shown 

for the two cases. 



+ 1 - 2/A - 



-,u(u 

J 



Ff^^ ± XAVA [ 1 - 



(11) 



The expressions for F2{x) and F^[x) are obtained as: 



F^{XA,Q^)=2 / d'r 



d?p M 



dp^ 5;^^"*™"(/, p, A:f,„)F2"^"*™"(a;jv, Q 



dp' 5r''*°"(p°, P, A:^,p)Fr*°"(^^, Q') 

2xnpI 



^ -1 + 



(12) 



i^3^(^A,Q')=2 / d'r 



d^p M 
(27r)3E(p) 



rf/ 5r*°"(p°, P, A;^,,)F|^''*°'^(x^, Q' 



J) cneutron^^O „ ^ ipneutron ^ Q^) 



+ / rfp'^5™"(/,p,fc^,„)F3^ 



P°7 - 
(pO _ p^^)^ 



(13) 



where 



7=^=1 + 



,2; AT 



(14) 



III. RESULTS AND DISCUSSION 



Using Eqs. (fT2|) and (fT3|) . we have calculated the and F^ structure functions for the 
lead nucleus with target mass correction (TMC) [28| and CTEQ6.6 parton distribution 
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FIG. 4. (Color online) Ratio R(x,Q2)=_£||^ with full calculation for = 5 GeV^ and = 50 
GeV^ using our base equation with TMC (dashed line) and the full model(solid line) at NLO. 
Here the results for the ratios R(x,Q'^)= ^pw^, losF^ ' 208^3^ 208^3^-= ^^^^ shown 
for the two cases. 



functions (PDFs) at the Leading-Order (LO) jSlJ. We call this as our base result 



hereafter 



we include pion and rho cloud contributions in F.f following the model of Ref. 26| and the 
shadowing corrections in and 13|, l29|, which we call as our full calculation (Total). 
These prescriptions have also been discussed in brief in our earlier works 15|, |25| . In Figs. ([1]) 
and (12]), we have presented the results for the ratios Ri{x, Q ) = 2osf^ 208]^^ ' " ~ ^'"^^ 
at = 5 GeV^ and compared them with the results obtained from the phenomenological 



studies of Hirai et al. 



161 and Eskola et al. 



20|. The details of our calculations in C 



and ^^Fe have been given in Ref. ISj- Here ^^C is a pure isoscalar nucleus while ^^Fe is a 
nonisoscalar one. The results for the ratios may be measured by the ongoing MINERz/A 
experiment as they are planning to perform the cross section measurements in a wide range 
of neutrino energies using ^^C, ^^Fe and ^ospb nuclear targets. We find that the results for 
the ratios obtained by using our base calculation and full calculation are different at low 
x<0.5 which is due to the fact that this is the region where the mesons cloud contribution 
dominates while for x>0.5, the only nuclear effects are Fermi motion and binding energy 
corrections. Furthermore, we find that the nature of the ratios for F2 and F3 are not the 
same while in the phenomenological studies they are almost the same in the present 

studied region of x. We also observe that the EMC effect is more prominent for ratios of 




Q^GeV^) 



Q^GeV^) 



FIG. 5. (Color online) F2{x,Q'^) vs in ^'^^Pb calculated using our base equation with TMC at 
LO, the full model at LO and NLO. 

structure functions between a heavy nucleus such as ^osp]^ g^j^^ ^ much hghter one hke ^^C 
than for ratios between a heavy nucleus (^'^'^Pb) and a medium-size one like ^^Fe. This fact 
can be observed clearly in Figs. ([3]) and @. 

In Figs.® and @ we have presented the results for ^po, ^pm, ^5^' ri^^ 



(z = 2, 3) at Q2 = 5 GeV^ and = 50 GeV^ using CTEQ PDFs at NLO [31]. These 
results are presented for base as well as full calculations. The deuteron structure functions 

n 

have been calculated using the same formulae as for ^^C in Ref. |15l|. but performing the 
convolution with the deuteron wave function squared instead of the spectral function. We 
have used the parametrization given in Ref. 



for the deuteron wave function of the Paris 



N-N potential. 
The results 

Figs. ([5]) and ([6]) where we have also shown the CHORUS experimental data [6|, l7|] for a wide 



The results for and structure functions in the lead nucleus have been shown in 
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FIG. 6. (Color online) xF3{x,Q'^) vs in ^ospb calculated using our base equation with TMC at 
LO, the full model at LO and NLO. 

range of x and Q^. The effect of shadowing is about 5-7% at x=0.1, Q^=l-5 GeV^ and 1-2% 
at x=0.2, Q^=l-5 GeV^ which dies out with the increase in x and Q^. In the case of there 
are pion and rho cloud contributions. The pion contribution is very dominant in comparison 
to the rho contribution. Pion contribution is significant in the region of 0.1 < x < 0.4. 
Thus, we find that the shadowing correction seem to be negligible as compared to the other 
nuclear effects. It is the meson cloud contribution which is dominant at low and intermediate 
X for F2. In these figures we also show the results of our full calculation at NLO. We can 
observe clearly that for medium Bjorken x ^ 0.6 there is no difference between the base 
calculation and the full model at LO which is because in this x-region the mesons cloud 
contribution has a negligible effect on the structure functions. We find an overall better 
agreement between the NLO calculation and the experimental data. 

We must point out here that in the case of lead nucleus the correction due to nonisoscalar- 
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FIG. 7. (Color online) -g^j^ vs y at different x for z^^ induced reaction in Pb at 



25 



GeV using CTEQ 



3l[ PDF at LO. Isoscalar full model: dashed line; nonisoscalar one: solid 



line. Full model without pion, rho and shadowing with isoscalarity: dotted line; nonisoscalarity: 
dashed-dotted line. 



ity is large and therefore the differential cross sections obtained by using the prescription 
given in our earlier paper 15| for the isoscalar case is not valid here. To observe the effect of 
nonisoscalarity, we have calculated -g^^ given in Eg. (HID using Eqs. (ll2p and ( !T3|) for ^osp]-, 
at iJj,^ p^=25 GeV and using Eqs.(13) and (14) of Ref. [15] treating lead to be isoscalar, and 
the results for and induced processes are respectively shown in Figs. ([7]) and (E]) for the 
full model as well as with the base term only. The cross sections have been calculated using 



CTEQ 3l| PDFs at LO. The difference is small at low x but increases with the increase in 
the values of x. For example, we find that the nonisoscalar correction to be around 4-5% at 
x=0.175 which becomes 10-12% at x=0.65. 

In Figs. (ED and we have shown the results for in ^ospb at E^^=25 GeV and 

70 GeV respectively. The calculations for the the double differential cross sections have 
been performed for > 1 GeV^. These results are presented for our base and full 
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FIG. 8. (Color onine) -^-^^ vs y at different x for induced reaction in ^'^^Pb at Ep„ = 25 
GeV using CTEQ 



311 ] PDF at LO. Isoscalar full model: dashed line; nonisoscalar one: solid 
line. Full model without pion, rho and shadowing with isoscalarity: dotted line; nonisoscalarity: 
dashed-dotted line. 



calculations at LO and the full calculations at NLO. Similarly in Figs. (ITT] ) and (fT2!) . we 

induced by antineutrinos in 208p]^ ^-^ _25 GeV and 70 



have shown the results for 



1 



E dxdy 



GeV respectively. It is difficult to disentangle the main effects that make the curves being 
different because in the cross section both structure functions are mixed. However, we can 
say that for neutrinos there is a clear difference between LO and NLO calculation while for 
antineutrinos this difference seems to be smaller. We also find that the results improve (in 
comparison to the experimental data) when we add the contribution of mesonic degrees of 
freedom to the base calculation and even more when we perform the calculation at NLO. 



IV. CONCLUSION 



To conclude, in this work we have studied nuclear-medium effects in the structure func- 
tions F^{x,Q^) and F^{x,Q'^) in lead nucleus using a many-body theory to describe the 
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FIG. 9. (Color online) -g^^ vs y at different x for u^{Eu^ = 25 GeV) induced reaction in ^°^Pb. 
Dotted line is the results using Eq(TT]with TMC. Solid (Dashed line) line is full calculation at NLO 
(LO). The experimental points are from CHORUS [7]. 



spectral function of the nucleon in the nuclear- medium for all . The local density ap- 
proximation has been used to apply the results for the finite nuclei. The use of th e sp ectral 



function is to incorporate Fermi motion and binding effects. We have used CTEQ [3l|] PDFs 
in the numerical evaluation. Target mass correction has been considered. We have taken 
the effects of mesonic degrees of freedom, shadowing, anti-shadowing in the calculation of 
F2 and shadowing and anti-shadowing effects in the calculation of F^. We have found that 
the mesonic cloud (mainly pion) gives an important contribution to the cross section. These 
numerical results have been compared with the experimental observations of CHORUS 
Collaboration. Using these structure functions we obtained differential scattering cross sec- 
tion for lead and compared the results with the CHORUS j?] observation. The results for the 
structure functions in carbon and iron {l5|, helium and oxygen have also been used to study 
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FIG. 10. (Color online) vs y at different x for i^ii{Eu^ = 70 GeV) induced reaction in ^°^Pb. 

Lines and points have the same meaning as in Figj9j 



the ratio ^ir^ (i=2,3), where for deuteron we have used 



the same formulae as for ^^C in Ref. [l5|, but performing the convolution with the deuteron 
wave function squared instead of the spectral function. These results have been compared 
with some of the phenomenological studies by Hirai et al 16| and Eskola et al. {2^. We 



find that the effect of the nuclear-medium is also quite important even for the deep inelastic 
scattering. Furthermore, we observe that the ratio of the structure functions for the different 

n 

nuclei are not the same. MINERz^A [Hi would be able to measure these ratios as they are 
using deuterium, helium, carbon, iron and lead targets and our results may be quite useful 
in their analysis. 
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FIG. 11. (Color online) -^^^^ vs y at different x for v^iEy^ = 25 GeV) induced reaction in ^°^Pb. 
Lines and points have the same meaning as in Figj9j 
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FIG. 12. (Color online) vs y at different x for Uf,{Eu^ = 70 GeV) induced reaction in ^ospb. 

Lines and points have the same meaning as in Figj9j 
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